Lateral diffusion of a protein on a fluctuating membrane 
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Measurements of lateral diffusion of proteins in a membrane typically assume that the movement of the pro- 
tein occurs in a flat plane. Real membranes, however, are subject to thermal fluctuations, leading to movement 
of an inclusion into the third dimension. We calculate the magnitude of this effect by projecting real three- 
dimensional diffusion onto an effective one on a flat plane. We consider both a protein that is free to diffuse in 
the membrane and one that also couples to the local curvature. For a freely diffusing inclusion the measured 
projected diffusion constant is up to 15% smaller than the actual value. Coupling to the curvature enhances 
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diffusion significantly up to a factor of two. 
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Introduction 

o . 

, During the last two decades significant progress has been made regarding novel microscopy techniques, that are mainly 
based on the observation of single or few molecules. With the aid of single particle tracking, photonic force microscopy, or 
fluorescence-based single-molecule microscopy one is now capable of studying the dynamic properties of lipids or proteins in 
membranes with positional accuracies smaller than 40 nm and a time resolution that can be as small as tens of microseconds, see 
the review 1 1] and references therein. These techniques together with other methods like fluorescence correlation spectroscopy, 
reviewed in |2J], have made it possible to measure lateral diffusion of single molecules in a membrane very accurately giving 
insight into the organisation of biomembranes. Although these high resolution methods have become standard practice, it is 
O [ typically overlooked that a membrane is soft and therefore subject to thermal shape undulations. Diffusion coefficients extracted 
from experimental data usually correspond to projected diffusion in a flat plane. The fluctuations of the membrane, however, lead 
t—{ \ to a three-dimensional motion of the diffusing protein, but unlike free three-dimensional diffusion the particle is confined to the 
i membrane. Since the fluctuation spectrum of the membrane can be affected by external parameters, like e.g. temperature, osmotic 
pressure or pH differences 121^!, these high accuracy methods should make the apparent change in diffusion experimentally 
measurable. 

In this paper, we are interested in how big the difference between the actual intramembrane and the measured projected 
diffusion constant is. We study both the case of a protein, which is free to diffuse 0101, and a curvature-coupled protein with (or 
ly-j ■ inducing) a spontaneous curvature. Most studies of membrane-inclusions within various models, which include rigid inclusions 
of various shapes, or proteins with a large domain outside the membrane, analyse the static interaction with the membrane and/or 
the interaction between inclusions , and have not yet addressed the dynamical issue of diffusion. 

a : 

I Free diffusion 

First, we regard the simplest case of a protein free to diffuse in the membrane. We assume that the diffusion is not influenced 
q , by the local shape of the membrane, therefore apart from the inclusion being confined to the membrane there are no interactions 
between membrane and protein. To describe diffusion on a curved surface, the Laplace operator in the diffusion equation for 
. a plane needs to be replaced by the Laplace-Beltrami operator. If the position of the surface f is expressed in the Monge 

■ representation, i.e. f = (x, y, h(x, yj), the resulting Smolouchovski equation is 01 

^ : 

dP(x,y,t) 1 f, , 9 ,d 2 P , l7 ,d 2 P , , d 2 P 
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D is the diffusion constant on the curved 



with h x = dh(x, y)/dx -for the other subscripts accordingly- and g = 1 + h 2 

surface, P(x, y, t)dxdy is the probability to find the diffusing particle in the area element dxdy at point f . The membrane is 
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subject to thermal fluctuations. Similarly to the calculation of diffusion within the Zimm model fl3ll we introduce a preaveraging 
approximation: Instead of using the prefactors, that contain h x , h y , h xx , h yy , h xy , and therefore explicitly depend upon position 
and time, we replace them by their thermal averages (...). This approximation is valid if the time scale of protein diffusion is 
much larger than that of membrane shape fluctuations. Due to their asymmetry averages like (h x h y /g) and many others vanish. 
Only (h 2 x /g) = (h 2 /g) are non-zero contributions. This considerably simplifies the diffusion equation 



dP(x,y,t) 
dt 




(2) 



Average quantities of an isotropic membrane cannot be different for the x- or y-direction, and the effective diffusion constant 
£> pro j, that would be measured in the x-y-plane is rescaled to 
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(3) 



This relation has been derived previously by a slightly different approach |6] but not evaluated. 

We evaluate the rescaling factor D pm j/D using the classical Helfrich Hamiltonian Ho as a model for the membrane. It takes 
the following approximate form in the Monge representation (r = (x, y), V = (d/dx, d/dy)) 



H Q [h(r,t)}= ^d 2 r~(V a fc) 2 + 5(Vfc) 3 , 



(4) 



where k is the bending rigidity of the membrane and a an effective tension (3- To calculate the expression (1/g) = 
/ °°da (exp[— ag]) we introduce the functional integral 



Z{p) = V[h] exp 



d 2 rUv 2 h) 2 
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(5) 



where L is the length of the system in the x- and y-direction and (3 = 1/ (fc^T) the inverse temperature. The partition function 
Z of the membrane is given by Z(a = 0) and therefore (exp[— a(V/i) 2 ]) = /-Z(O). Using a Fourier expansion of spatially 
varying variables the calculation of i?(a) is straightforward and yields 



(exp[-a(V/i) 2 ]) = exp 
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The cutoff wave number q m ~ 1/a is given by the smallest length scale a present in the membrane. Because Lq m > 1 we 
replace the summation on the rhs by the integration over k. The ratio of projected and intramembrane diffusion constant may 
then be written as 
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This integral is a function of three dimensionless quantities Lq m , (3k, and (3a/ q^. An analytic solution is not found, but the fast 
decay of the integrand facilitates the numerical evaluation. 

For a numerical calculation of the scaling factor of the diffusion constant, we first need to analyse the size of the parameters 
(3, k, a, L, q m which go into eq. 0. Typical values of (3k for lipid bilayer membranes lie in the range between 5 and 50. The 
smallest length scale a of the system, that gives q m ~ 1/ a, is roughly the size of a lipid which is on the order of nanometres. The 
area L 2 of membranes studied in experiments can vary strongly from a few /im 2 to approximately 100 000 /im 2 . For iLq m ) 2 
we, therefore, regard the range from 10 6 to 10 11 . Typical values for the effective tension of a fluctuating membrane at room 
temperature are 10~ 6 mJ/m 2 < a < 10~ 3 mJ/m 2 , while rupture occurs for tensions on the order of a ~ 1 mJ/m 2 11411 . The 
whole range we reg ardis 10~ 7 s$ (3a/q 2 m ^ 10" 1 . 

First, we analyse the case of vanishing effective tension a = 0: In fig.^we show D pro j/D as a function of iLq m ) 2 for the 
given values of (3k. Overall we see that in the parameter ranges, which correspond to experimental conditions, the projected 




FIG. 1: Rescaling factor D pm j/D as a function of (Lq m ) 2 for FIG. 2: D„ mi /D as a function of rigidity (3k for the given values 

the different given values of (3k and a = 0. The softer the of (Lq m ) 2 . 
membrane, the stronger is the rescaling of the diffusion constant. 
The larger the regarded membrane, the stronger is the shift in D. 



diffusion constant is reduced by up to ~ 15%. As was to be expected we also find that a more rigid membrane leads to weaker 
rescaling of the diffusion constant and that an increase in membrane size, which leads to stronger fluctuations, causes the scaling 
factor to decrease. In fig.|2]we display D pro j/D as a function of rigidity (3k for the given (Lq m ) 2 . We find that an increase in (3k 
from 5 to 50 increases D pm j/D monotonically by approximately 0.1. For larger membranes, when the role of the fluctuations is 
more pronounced, D pm] /D is smaller and the rise with (3k stronger. 

We now consider the influence of an effective tension. The larger a the more expensive it is for the membrane to fluctuate, 
i.e. to build up gradients |Vft|. As a consequence, an increase in a will lead to a weaker reduction of D pm JD. This can be 
seen in fig.|5] where we display D pm j/D as a function of (3k and (3a /q 2 ^ for (Lq m ) 2 — 10 7 . For the smallest regarded effective 
tension and bending rigidity the reduction of the diffusion constant is the strongest with ~10%. An increase in (3k and (3a/q 2 n 
leads to a decreased difference between the projected and actual diffusion constant. For larger effective tensions and smaller 
rigidity the lines of constant D pm j/D appear almost linear. This behaviour can be extracted analytically from eq. Q in the limit 
of large systems 

D pm] /D* {l+ [I + + pK/(Pa/q 2 m )) /(Air (3k)]- 1 } /2, (8) 
which may be a more convenient expression for future reference than the full expression 0. 



Curvature-coupled diffusion 

In this section we regard a protein that interacts with the fluctuating membrane. The exact mechanisms, how a protein couples 
to a bilayer membrane are not understood, however there are a large number of possibilities that have been explored in theoretical 
calculations and simulations f7l l9l fiol [HI O . In this study we model a protein of radius a p at position R = (X, Y) in the 
projected plane to couple to a local curvature C v of the membrane via the elastic coefficient m. This situation could be caused 
either by the shape of the protein alone, or by a protein with a large extramembrane domain, that induces a spontaneous curvature 
in the membrane H2I1 . The energy Hi of this interaction, 

HrlHR)] = j (Vft/i(R) - C p ) 2 ~™ 2 p (V^(R)) 2 , (9) 

needs to be added to the membrane energy Hq of eq. @ to give the full energy H = Ho + Hi of the system. The second term 
of eq. accounts for the fact that there is no membrane where the protein is. If we are interested in the dynamics of the system 
we need to take the equations of motion both for the membrane and the particle into account. The dynamics of the membrane is 
expressed by 

dh(r)/dt = - f d 2 r'A(r, r') (SH /Sh(r') + 5Hi/6h(r')) + £(r) (10) 



with the kinetic or Onsager coefficient A(r, r') and random fluctuations £, which obey the fluctuation dissipation theorem. In 
Fourier space, the kinetic coefficient for a free, on average planar, membrane embedded in infinite space is given by A(k) = 
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FIG. 3: Factor D plo j/D as a function of /3k and /3a /q^ for 
(q m L) 2 = 10 7 . The reduction in diffusion constant is the weak- 
est for stiff membranes at high effective tension. 



FIG. 4: Ratio D cc /D of the diffusion constant of a curvature- 
coupled protein and intramembrane diffusion as a function of 
rigidity /3k and effective tension /3a /q^ for (Lq m ) 2 — 10 r . 



(Ar]k) 1 with the viscosity 77 of the fluid surrounding the membrane 1 14]. The dynamics of the inclusion follows the Langevin 
equation 

dR/dt = -Mproj V R W + C = Mproj [mC p + (mra 2 p - m)V|fc(R)] V R (Vf^R)) + C, (11) 
where £ is a random force acting on the inclusion with 

(C(*)> =0 and (Ci(*)Ci(*')> =2D Jwj M(t-f), (12) 

while /i pro j is a mobility, that is connected to the diffusion coefficient via the Einstein relation fcsT/i pl0 j = D pro j. Because we 
are only regarding the movement of the particle in the projected plane we need to use the projected diffusion constant £> pl0 j, 
which we previously related to the intramembrane diffusion constant D. The solution of the coupled eqs. JlQi and (II It defines 
the dynamics of the system. If we regard eq. (II Oi . it is clear that the solution will have two additive parts: the first is the result 
ho(r, t) of the membrane without inclusion, while the second hi (r, i) is a correction caused by the protein-membrane interaction 
energy Tit- To calculate the full dynamics of the particle ho + hi is plugged into eq. il It . If, however, the interaction energy 
is much smaller than the unperturbed membrane energy, Hi <C Tto, we may use the first order approximation of the Langevin 
equation (II It . which corresponds to using only the dynamics ho of an unperturbed membrane to calculate the protein dynamics. 

This approximate form of eq. il It is now used to calculate the diffusion constant of the curvature-coupled protein inclusion 
defined as D cc = lim t _ +00 (AR 2 )/4i. AR 2 is the projected squared distance in the X-Y-plane the particle has moved during 
time t. In the following we will drop the subscript of h and assume that it is the result of the unperturbed membrane. Using 
eq. <|TT) we write (R(f = 0) = 0) 
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drdrV p 2 roj {m 2 C p 2 ([V(V 2 ^t))][V(V 2 /*(t'))]) + 0(h 4 ). (13) 



For the last line we used eq. H2\ and remembered that averages of uneven powers of h vanish. Because contributions on the 
order of 0(h 4 ) are negligible, they are omitted in the following. To calculate the diffusion coefficient, we need the thermal 
average (M) = ([V(V 2 /i(t))][V(V 2 /i(0))]). We introduce the function Z (a) 



Z(a)= Jv[h]exp -Pj^yi Jd 2 k^^-^k 6 exp[-A(k)E(k)t]jh(k)h*(k) 



(14) 



with E(k) = nk 4 + ak 2 . The desired quantity is then determined by (M) = -gigj . Inserting the evaluation of eq. dl4> 
into eq. (II 31 and performing the two time integrals yields 



/*t,2, s\ .r^ 2 2^1 r m ,, 2 ^ 7 f exp\-A(k)E(k)t] - 1, 
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from which we derive the diffusion constant D cc for curvature-coupling as 



Mpro 



2 r 2 1 1m 
* 7T K z 



2(l + f^) 



■ arctan 



(16) 



Contrary to the last section, where we could give quantitative predictions for the relation between projected and free intramem- 
brane diffusion, this expression is subject to more model parameters and is less universal. However, it is still instructive to make 
a semi-quantitative estimate to gain insight into the influence of curvature-coupled protein diffusion. 

The factor [. . .] on the right of eq. (I16> . that contains only yj n/aq^, is always very close to one. Only towards the low end 
of the regarded range 10 ^ yj K/aq^ ^ 10 4 is this factor slightly reduced. Therefore, the prefactor of [...], which contains the 
elastic constant m, the spontaneous curvature C p , the rigidity k, the cutoff wave number q m , the viscosity 77, and the mobility 
/ip ro j, determines the magnitude of the diffusion. Because this prefactor is obviously positive, the diffusion of a protein is 
generally enhanced for the considered coupling to local curvature. 

The characteristic scale of the elastic constant m is the bending rigidity k of the membrane times the area na^ occupied by 
the protein. Thus we can write m = c\ n na^ with a constant c\ on order of unity. We choose c\ = 2. An axisymmetric 
inclusion, whose area on one side of a membrane with width d is twice that of the other side, produces a spontaneous curvature 



of C p = 4 /3d 1 14]. A good estimate for the cutoff wave number is q n 



la 



We further assume the radius of the protein to be 



a, p = 8a, and the width of the membrane d = 5a. These estimates ensure that the energy of the protein may still be regarded as 
a perturbation. Unlike the initial problem, where we could give the ratio D pmi /D, we now also need an estimate for the mobility 
Mproj = PD vm y Saffman and Delbriick 1 15] derived this expression for the mobility within a flat two-dimensional liquid layer 
bound by a surrounding three-dimensional fluid with viscosity 77 



// = [\n{yd/{r)a p )) - 7] /(4W). 



(17) 



In this equation, v is the viscosity of the liquid layer, i.e. the membrane, and 7 ~ 0.577 is the Euler number. This mobility 
is valid for vd/r)a p >> 1, which is the case for a typical lipid membrane. We use eq. (I17t with the previously calculated 
rescaling, /i pl - j = [i' D pm J D, for the estimate of the diffusion constant D cc . For a lipid membrane, the viscosity is on the order 
of v ~ lerg sec/cm 3 and the viscosity of water is r\ ~ 10~ 2 erg sec/cm 3 . This gives r\jv ~ 0.01. If we neglect the effective 
tension and use our estimates in eq. fl!6l >. we find: 



D cc (a = Q)/D ~ D pmi /D {1 + 0.81 D pmi /D} 



(18) 



For a typical ratio of D pm j/D — 0.95 we see that the diffusion of the protein is increased by ^68%. This estimate reveals that a 
coupling to local curvature can lead to a significantly enhanced diffusion coefficient. 

The ratio D cc /D for a > using D pm j/D from eq. @ is shown in fig.|4]as a function of (3a /q^ and (3k for (Lq m ) 2 = 10 7 . 
Surprisingly the weakest influence on the diffusion constant is found for small rigidity (3k and small effective tension (3a/ q^, 
i.e. for membranes with strong fluctuations. This follows from a partial compensation of the two effects discussed in this paper: 
while the curvature-coupled interaction between membrane and protein enhances diffusion, the fluctuations reduce the projected 
diffusion in the x-y-plane, cf. fig.|3] Starting from small rigidity and tension we see that the diffusion coefficient increases when 
the rigidity and effective tension are raised. Increasing (3a /q^ and (3k leads to weaker fluctuations in the membrane, which 
finally becomes unfavourable for the protein-membrane interaction. This becomes especially visible, in the case of constant (3k 
and increasing (3a/qf n : initially, an increase in the diffusion constant is observed followed by a subsequent decrease. 



Concluding perspective 

In summary, we have found that thermal fluctuations of biomembranes have considerable influence on the diffusion of proteins. 
When the protein is free to diffuse within the membrane, the projected diffusion constant, which corresponds to the quantity 
typically measured, is up to 15% smaller than that of the true intramembrane diffusion. The lower the bending rigidity or the 
effective tension and the larger the membrane, the stronger is this effect. Experimentally, it could be studied by observing the 
change in the projected diffusion constant upon osmotically induced swelling of an initially flaccid vesicle, which increases 
the effective tension 11411 . To gain insight into the influence of protein-membrane interactions, we studied the diffusion of 
an inclusion with a spontaneous curvature and found that this enhances diffusion significantly. When the projected diffusion 
constant is measured this leads to an interesting interplay: on the one hand fluctuations can lead to enhanced intramembrane 
diffusion, while on the other hand stronger fluctuations lead to a smaller projected diffusion constant. It is therefore possible 
that the measurement of the diffusion coefficient at constant rigidity and increasing effective tension reveals a maximum for a 
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certain effective tension. Future experimental studies should analyse the effect of changing membrane fluctuations on the lateral 
diffusion of proteins, since they could shed light on the local coupling mechanisms between proteins and lipids. 
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